Time-resolved photoemission by attosecond streaking: extraction of 
time information 
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Abstract. Attosecond streaking of atomic photoemission holds the promise to provide unprecedented information on the release time of 
the photoelectron. We show that attosecond streaking phase shifts indeed contain timing (or spectral phase) information associated with the 
Eisenbud-Wigner-Smith time delay matrix of quantum scattering. However, this is only accessible if the influence of the streaking infrared 
(IR) field on the emission process is properly accounted for. The IR probe field can strongly modify the observed streaking phase shift. 
We show that the part of the phase shift ("time shift") due to the interaction between the outgoing electron and the combined Coulomb 
and IR laser fields can be described classically. By contrast, the strong initial-state dependence of the streaking phase shift is only revealed 
through the solution of the time-dependent Schrodinger equation in its full dimensionality. We find a time delay between the hydrogenic 
2s and 2p initial states in He+ exceeding 20 as for a wide range of IR intensities and XUV energies. 

PACS numbers: 32.80.Fb, 42.50.Hz, 42.65.Re, 31.15.A- 



Attosecond streaking is one of the most spectacular 
applications within the emerging field of attoscience [1- 
7]. It has developed into a powerful tool to achieve time- 
resolution on the sub-100 attosecond time scale [8-11]. 
Streaking is based on a variant of a pump-probe setting 
with an extreme ultraviolet (XUV) pulse of a few hundred 
attoseconds serving as pump and a phase-controlled few- 
cycle infrared (IR) pulse as probe. Electrons emitted in 
the presence of an IR field are accelerated to different final 
momenta and energies depending on the value of the vector 
potential at the release time. Thus, time information is 
mapped onto the energy axis in analogy to conventional 
streaking. First proof-of-principle implementations were 
the direct time-domain measurement of the life time 
of the Xe(4p _1 ) hole by Auger decay of ~ 8fs [8] 
complementing equivalent spectral information, as well as 
the measurement of the time-dependent electric field of 
the IR light wave [9]. More recent and more advanced 
applications addressed the delayed photoemission from 
core levels relative to conduction band states of a tungsten 
surface due to the increased pathlength and inelastic 
scattering processes [10]. Very recently, a time delay of 
atomic photoemission from the 2s shell relative to the 2p 
shell of Neon as small as 21 attoseconds, i.e. less than one 
atomic time unit has been extracted [11]. This time scale 
is one order of magnitude shorter than the XUV "pump" 
pulse and two orders shorter than the oscillation period T 
of the IR probe pulse. First theoretical model calculations 
[11, 12] could qualitatively account for a delay of the 2p 
relative to the 2s electron, falling, however, short by more 
than a factor of 3 to quantitatively reproduce the ~ 20 as 
delay. 



As the resolution of the phase of the streaking field 
reaches such a small fraction of an optical cycle that 
the associated time resolution approaches the atomic 
time scale (a classical electron in the ground state of H 
travels 0.5 A during ~ 24.2 as), fundamental conceptual 
questions arise as to the physical information to be 
accessed. More generally, the question is posed which novel 
information can be extracted in the time domain which 
cannot be gained by high-resolution measurements in the 
spectral domain. For complex systems in the presence 
of incoherent processes (such as electron transport near 
surfaces [13]), time domain studies of relaxation and 
transport can obviously provide novel information on 
competing decoherent pathways not easily accessible in the 
spectral domain. For fully coherent quantum processes 
where time and spectral domain descriptions are, in 
principle, linked by a Fourier transform, the extraction 
of novel information is of considerable interest. 

In this paper we show that attosecond streaking opens 
the pathway towards a direct observation of matrix ele- 
ments of the Eisenbud-Wigner-Smith (EWS) time delay 
operator [14, 15] down to a remarkable attosecond level. 
The EWS time delay is a measure for the spectral variation 
of the scattering phase. For later reference we note that 
the EWS time delay pertains to scattering in short-range 
potentials. Long-range Coulomb interactions are beyond 
the realm of standard EWS theory and require special care. 
As we will show, the extraction of the EWS time delays 
from streaking phases is highly non-trivial as the IR probe 
field and the long-range Coulomb field distort the spec- 
tral phase information to be extracted. The distortion can 
occur both in the entrance channel of the atomic bound 



state to be photoionized as well as in the exit channel 
of the continuum electronic wave packet leaving the ionic 
core. Remarkably, the latter can be accurately accounted 
for within a single-particle description by classical dynam- 
ics. A classical-trajectory Monte Carlo (CTMC) simula- 
tion is found to be in excellent agreement with the nu- 
merical solution of the time-dependent Schrodinger equa- 
tion (TDSE). For the entrance channel distortion we find 
a strong state and angular momentum dependence. In a 
proof-of-principle simulation we show that, after account- 
ing for IR field distortion in the entrance and exit channels, 
the sensitivity of the streaking phase to short-ranged ad- 
mixtures can reach the level of a few attoseconds. Atomic 
units will be used throughout the text, unless indicated 
otherwise. 

We numerically solve the full three-dimensional TDSE 
within a one-electron (or independent particle) framework, 
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where Fxvv(t) is the (linearly polarized) electric field 
of the attosecond pump pulse (huixuv ~ 100 eV), 
the intensity of which is sufficiently weak (ixuv ~ 
10 13 W/cm 2 ) as to allow for a perturbative treatment, 
and -Fnt(i) is the electric field of the few-cycle IR probe 
pulse (fiwiR 1.5 eV) with sufficient intensity in the non- 
perturbative regime such that the energy modulation by 
the streaking field provides energy (and time) resolution 
but remains below the level where tunnel ionization 
becomes prevalent (7ir « 10 10 — 10 12 W/cm 2 ). The electric 
fields are related to the corresponding vector potentials by 
F(t) = —M^A{t). For the spherically symmetric atomic 
binding potential V(r) we employ pure Coulomb potentials 
Vc(r) = — — , as well as short-ranged model potentials. 

Our computational method for solving Eq. 1 is based 
on the well-established pseudo-spectral split-operator 
method as described in [16]. The time-evolution is 
carried out using a split-operator algorithm where the 
evolution in the Coulomb field is performed in the 
energy representation of the wave function and the action 
of the external field in length gauge is calculated in 
coordinate space. The energy eigenstates are obtained 
by diagonalizing the field free Hamiltonian in a finite 
element discrete variable representation basis (cf. [17- 
20]). In order to disentangle channel distortion effects 
by the probe field from the ionization dynamics driven 
by the XUV field, we also consider a constrained TDSE 
("restricted ionization model" [21, 22]) where the wave 
function ^f(f,t) is decomposed into the initial state 3>o(f) 



and the continuum part <f>(f, t) created by absorption of 
a single photon from the XUV pulse. The IR field is 
constrained to act only on t). 

In our classical simulations we propagate electrons in 
the combined Coulomb and laser fields using a standard 
variable step-size Runge-Kutta algorithm. 10 7 initial con- 
ditions were chosen either from the microcanonical ensem- 
ble for the binding energy of the level with principal quan- 
tum number n or from its (n, £, m) dependent subensem- 
bles [23] where I and m denote the orbital and magnetic 
quantum numbers, respectively. At the ionization time r a 
randomly chosen energy from the energy spectrum of the 
XUV pulse (Gaussian distribution with a FWHM of 4eV) 
is transferred to the electron by adding a momentum in 
random direction. The ionization probability is taken to 
be proportional to the XUV intensity. 

Starting point of the attosecond streaking technique 
[3] is the assumption that the asymptotic momentum of a 
photoelectron ionized by an XUV pulse of duration of a few 
hundred attoseconds is shifted by the instantaneous value 
of the vector potential Air(t) at the moment of ionization, 
i.e. arrival in the continuum, 



Pf( T ) =Po~ 4tr(t) 



(2) 



relative to the unperturbed asymptotic momentum of the 
photoelectron po = [2(wxuv — emitted by a photon 

with energy wxuv and initial binding energy £\ (in a.u.). 
For temporally well-defined IR electric streaking fields 
Fin(t) and vector potentials 



(3) 



Air(t) = - / F m (t)dt 



the momentum shift and, hence, energy shift of the 
photoelectron allows thus, in principle, a mapping of 
emission time r onto energy ("streaking") and therefore 
the determination of emission time or, more precisely, of 
the phase <f>s = wirt of a streaking spectrogram. An 
accuracy of the order of ten attoseconds or ~ 0.05 rad. 
can be reached for a streaking IR field with an oscillation 
period of the order of one femtosecond (T s» 2.6 fs for 
Air = 800 ran). The physical interpretation and analysis 
of streaking spectrograms and derived phases or times 
require, however, a detailed inquiry into the assumption 
underlying Eq. 2. 

The fundamental assumption of the streaking field as 
a probe is that the IR field does not distort the system to 
be probed. Since for easily resolvable energy shifts IR fields 
with intensities of the order 7ir « 10 11 — 10 12 W/cm 2 are 
needed, distortion effects generally cannot be neglected. 
Unlike for the weak XUV pump pulse, effects beyond 
lowest-order perturbation theory must be considered for 
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Figure 1. Streaking spectrograms for an 800 nm IR laser field with 
a duration of 3fs and an intensity of 10 12 W/cm 2 . The graphs show 
the final momentum distribution in "forward" direction of the laser 
polarization axis for H(ls) and He+(ls) initial states and an XUV 
photon energy of 80 eV. The solid white lines are the first moments 
of the electron spectra. The dashed white vertical lines indicate the 
shift of the central minimum relative to the vector potential. The 
colour scale is in arbitrary units. 

the IR field. As we will show in the following, the near 
independence of 4>s on Jtr within the accessible range does 
not contradict the notion that (j>s is strongly affected by 
the IR field. Moreover, the assumption of a "sudden" 
transition to the undisturbed momentum po caused by 
photoabsorption (Eq. 2) followed by the momentum shift 
by the IR field of the remaining streaking pulse requires 
scrutiny when we aim for the resolution of electronic 
dynamics on the attosecond scale. As the wave packet 
recedes from the ionic core, it propagates for ~ 100 as in 
the atomic (ionic) potential with a local momentum p(r) 
rather than the asymptotic momentum pq. Neglect of the 
deviation of the local from the asymptotic momentum is 
in line with the strong-field approximation underlying the 
original attosecond streaking model [3] but is no longer 
valid for streaking in the Coulomb field. 

Photoionization spectrograms for ionization of the 
H(ls) and He+(ls) states by an 80 eV XUV pulse with 
a duration of 200 as (FWHM of the Gaussian intensity 
envelope), streaked by an 800 nm IR laser field with a 
duration of 3fs (FWHM of the cos 2 envelope) and an 
intensity of 10 12 W/cm 2 (as well as the first moment of 
the final electron momentum distribution, white lines) 
indicate the level of phase (or time) resolution achievable 
(Fig. 1). By fitting the streaking curves to the analytic 
form of the IR vector potential Am(t + ts), we obtain 
an absolute time shift ts (with respect to the vector 
potential) of —6.9 as for the ionization from the H(ls) state 
and —37.1 as for the ionization from the He + (ls) state, 
yielding a relative delay of At$ = 30.2 as between the 
two emission channels. The sign convention for ts ensures 
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Figure 2. Temporal shifts t$ extracted from quantum mechanical 
streaking simulations (H(ls): red circles, He + (ls): blue circles), 
classical streaking simulations (H(ls): red squares, He + (ls): blue 
squares), and for comparison, the EWS time shift <ews = dip/dE 
applied to the Coulomb phase (H + : red solid line, He 2+ : blue dashed 
line) . 

that positive values correspond to delayed emission, i.e. the 
electron "feels" the vector potential at a later time. The 
spectrograms are obtained with a full numerical solution 
of the 3D TDSE without any approximations. A negative 
time shift, or time advance, ts < 0, relative to the 
vector potential can be observed over a wide range of 
final electron energies (Fig. 2) corresponding to varying the 
XUV pulse frequency w X uv- All results are obtained from 
spectrograms taken in "forward" direction with respect to 
the laser polarization axis and an opening angle of 10 deg. 
(in the case of s-like initial states there are, however, no 
significant changes for the "backward" spectrograms and 
different opening angles). The error bars for the retrieved 
time shifts correspond to the estimated standard error 
from the nonlinear least-squares fit of the analytic vector 
potential to the calculated spectrograms. 

The observed time shifts ts (or phase shifts </>s) are 
not directly related to the EWS-like time shift given by 
the energy derivative of the spectral phase ip, i.e. the group 
delay [24] 



also shown in Fig. 2. For the current case of 

photoionization from Is initial states in Coulomb 

potentials the phase ip would be given by the Coulomb 
phase 

(p = a e = arg[r(j? + 1 - iZ/k)] (5) 

where I = 1 is the angular momentum of the free electron, 
k its wavenumber, and Z is the charge of the remaining ion. 
Even the sign of £ews predicted by Eqs. 4 and 5 is different 
from the streaking time shifts ts (Fig. 2). Obviously, £ews 
is overshadowed by the distortion of the outgoing Coulomb 



continuum electron by the IR field. It should be noted 
that the EWS time shift is only well-defined for short- 
ranged potentials. It is therefore no surprise that Eq. 5 
is unrelated to what is observed in terms of the streaking 
phase. 

Since classical-quantum correspondence holds for an 
unbound electron in both the Coulomb field (Rutherford 
scattering) and a laser field (Volkov states), it is suggestive 
to apply classical dynamics also to the motion of the 
outgoing electron in the combined Coulomb and IR fields 
in a streaking setting. Accordingly, the final momentum 
of the electron on its trajectory taking off at r(r) s» near 
the ionic core is given by 



(6) 



p f ( T )=p + / dta F c {f(t)),F m (t) 



In the limit of vanishing Coulomb field, Fc(t) = 0, 
the acceleration a reduces to a[Fm,(t)} — and 
Eq. 6 becomes equal to Eq. 2. For vanishing streaking 
field, i<iR,(i) = 0, Eq. 6 describes the approach of the 
instantaneous local momentum p(r) to the asymptotic 
momentum p along the Kepler hyperbola as t — > oo (or 
r — > co). The Coulomb corrected mapping from emission 
time to momentum shift thus becomes 



Pf( T ) =Po~ A m (r)+ 

dt (a[P a (T%t)),P m (t) 



F m (t) 



(7) 



Eq. 7 represents the classical realization of Coulomb- 
laser coupling (cf. [25, 26]) in the exit channel through the 
modification of the trajectory probing the Coulomb field 
by the simultaneous presence of the IR field. It treats the 
IR field and the Coulomb field non-perturbatively and on 
equal footing. One can easily show that the correction 
term relative to the standard streaking relation Eq. 2 
yields, to leading order, a term ~ -Fir(t) linear in the 
streaking field but phase shifted relative to Air(t) by ir/2. 
Consequently, Eq. 6 becomes 



p f (r) KipQ-a- A m (r + t s ) , 



(8) 



with ts the Coulomb-laser coupling induced time shift. 
Remarkably, ts is, to leading order, independent of the 
amplitude of the IR field since it depends on the interplay 
between F m and A m (cf. Fig. 3). 

Implementation of Eq. 6 within our classical trajec- 
tory Monte Carlo (CTMC) simulation gives streaking time 
shifts in an almost perfect agreement with the full TDSE 
over a wide range of final-state energies for both H(ls) 
and He + (ls) initial states (Fig. 2) and over two orders of 
magnitude in intensity (Iir rs 10 10 — 10 12 W/cm 2 ) of the 



streaking field (not shown). The approximate indepen- 
dence of ts on the streaking field intensity following Eqs. 
6 and 8 may mask the origin of ts to be, indeed, due to the 
Coulomb-laser coupling. We note that the present classi- 
cal phase shift due to Coulomb-laser coupling appears to 
agree better with the full TDSE than the semi-classical 
estimates based on the eikonal approximation [27] which 
may be due to the perturbative approximation involved. 

A second class of distortion effects refers to the 
entrance channel, i.e., to phase shifts of the initial state 
by the IR field. Clearly, those phase shifts are expected to 
be state dependent and are outside the scope of classical 
dynamics. We restrict our investigation to low-lying states 
of He + (n£m) which are stable against tunnel ionization 
by the IR field (J m < 10 12 W/cm 2 ). Within the full 
TDSE calculation we find remarkably strong initial-state 
dependent streaking field distortions. Most notably, the 2p 
substates aligned along (m = 0) and perpendicular (m = 
1) to the polarization axis acquire a relative streaking 
phase shift of Ats ~ 15 as, and the 2s relative to the 2po of 
Ats ~ 10 — 20 as. This pronounced streaking phase shift 
as a function of angular momentum differs from recent 
calculations [27, 28] which did not find a shift for states 
with (initially) vanishing dipole moment. 

By performing constrained TDSE calculations ("re- 
stricted ionization model") in which the initial-state in- 
teraction with the IR field is suppressed but continuum 
states fully evolve under the influence of the IR field, we 
find that the relative phase (time) shifts between 2s, 2p , 
and 2pi, for the most part, disappear. Thus, a state- 
dependent exit channel distortion can be ruled out. We 
find that also the initial-state dependent phase (time) shift 
is, within the range of useful streaking field intensities 
(/ir sa 10 10 — 10 12 W/cm 2 ), remarkably insensitive to the 
value of /ir (Fig. 3) closely following the behaviour ob- 
served for the exit channel distortion. This apparent in- 
sensitivity to /ir should not mislead to the conclusion that 
distortion effects by the IR field on the atomic dynamics 
can be excluded. 

The strong variation of ts with angular momentum 
£ is not specific to the degeneracy of the hydrogenic 
Hc + (n = 2) manifold. By choosing an atomic model po- 
tential that breaks the SO (A) symmetry we have verified 
that the streaking time shift persists for non-hydrogenic £ 
manifolds, thereby excluding degeneracy effects as the ori- 
gin. Moreover, transient inter-shell coupling was found to 
dominate dynamical polarization over intra-shell mixing. 
One remarkable observation of the present results (Fig. 3) 
is that the relative time delay between an emitted 2p elec- 
tron with 80 eV final energy relative to a 2s electron with 
a final energy of 55 eV yields Ats ~ 20 as in striking agree- 
ment with the first experimental data for neon at the same 
final-state energies [11]. Clearly, this agreement may be, 
in fact, fortuitous as the electronic structure of neon is not 
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Figure 3. Streaking time shifts tg extracted from quantum 
mechanical streaking simulations for different initial states. The 
results are obtained from spectrograms taken in "forward" direction 
with respect to the laser polarization axis and an opening angle of 
lOdeg. . 

properly accounted for. To the extent, however, that the 
dominant contribution to Atg results from £ dependent 
entrance and quasi-classical exit channel IR distortions, 
the present results may, indeed, account for a considerable 
fraction of the delay observed. 

So far, the present results quantitatively demonstrate 
the profound influence of the IR field on the observed 
streaking time shift. In the following, we will show that 
atomic photoemission times can be, indeed, accessed by 
streaking, provided that IR distortion effects can be cither 
excluded or accounted for. We illustrate the power of 
attosecond streaking with the help of two prototypical 
examples. 

We first consider an isotropic Is- like initial state 
deeply bound at £j sa — 54.4eV by a Yukawa- type short- 
ranged model potential with screening radius a, 

W(r) = --e-». (9) 
r 

In this case, both initial-state dynamical polarization by 
the IR field and exit channel Coulomb-IR field coupling 
can safely be excluded. We calculate the EWS time 
delay or, equivalently, the energy variation of the spectral 
phase, £ews = dipy/dE, for a wide range of final state 
energies and different combinations of charge Z and atomic 
screening radii a such that £{ remains constant. We 
find near perfect agreement with the corresponding time 
shift ts extracted from streaking spectrograms down to 
~ 15 eV final electron energy (Fig. 4). Interestingly, 
£ews varies non-monotonically with a. Both time delays 
(^ews > 0) and time advances (^ews < 0) relative to 
the field ^4rR can occur. A physical realization of this 
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Figure 4. Temporal shifts tg extracted from quantum mechanical 
streaking simulations for the short-ranged Yukawa potential Vy with 
screening lengths a = 0.5 (squares), a = 1 (triangles), and a = 2 
(diamonds) a.u. . The streaking IR laser field has a wavelength of 
800 nm, a duration of 3fs, and an intensity of 10 12 W/cm 2 . Points: 
streaking shifts tg, lines: <ewS (Eisenbud-Wigner-Smith "delays"). 

scenario would be photodetachment of negative ions. This 
is, however, complicated by the fact that typical streaking 
field intensities would rapidly destroy the weakly bound 
negative-ion state. 

As a second scenario we consider a short-ranged 
admixture to the Coulomb potential. In this case, the 
exit-channel distortion by the IR field is present but can 
be accurately accounted for (Eq. 7). For a short-ranged 
admixture of the form 

(r— r ) 2 

SV(r) = Vq • e 5^ (10) 

to the tail of the Coulomb potential the residual EWS 
time shift can be determined. 5V(r) is chosen such that 
the initial state wave function and binding energy re- 
main undistorted. This requires that ro > (r)i s and 
|(5V(r)| <C |Vc(r)| such that the photoelectron energy 
distribution remains undisturbed. We compare now the 
EWS time shift (averaged over the spectral width of the 
XUV pulse) induced by Eq. 10 with the observed differ- 
ence in the streaking time shift, At$ — ts — (ts)c> after 
the Coulomb-laser coupling induced shift (ts)c has been 
subtracted. We find remarkable agreement for different 
locations ro of the local perturbation and both signs of 
Vq (Fig. 5). Even more remarkable is that the emission 
time shifts on the single-digit attosecond scale become ac- 
cessible. A possible physical realization of this scenario 
would be photoemission from the central atom of endohe- 
dral Ceo where the outgoing electron traverses the short- 
ranged shell potential of the cage. 

In summary, we have identified on the one-electron (or 
independent particle) level considerable state-dependent 
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Figure 5. Streaking time shifts Atg relative to the pure Coulomb 
case with Z = 2 for a short-ranged admixture Eq. 10 (Vrj = ±0.1 a.u., 
d = 0.33 a.u.) at different positions tq. The points are the time 
shifts as extracted from full streaking calculations for an XUV photon 
energy of 90 eV, an 800 nm IR laser field with a duration of 3 fs, and 
an intensity of 10 12 W/cm 2 , the lines are tBWS f° r the short-range 
admixture to the Coulomb potential, i.e. the Eisenbud-Wigner-Smith 
time "delay" . The results of the EWS time delay were convoluted 
with the pure XUV spectrum to account for the spectral spread of 
the electronic wave packet. 



time shifts that can be observed in attosecond streaking 
and which are of quantum mechanical origin. In addition, 
we have identified large time shifts which result from the 
coupling between the IR streaking field and the Coulomb 
field which depend on the final energy of the free electron 
and can be accounted for classically. The Eisenbud- 
Wigner-Smith (EWS) time shift (or energy variation of 
the scattering phase) is found to be accessible by streaking 
only if both initial-state dependent entrance channel and 
final-state exit channel distortions are properly accounted 
for. For such a scenario we have shown that time delays 
on the single-digit attosecond scale due to short-ranged 
potentials are in reach. 
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